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Abstract 

Neutrino emission due to the pair breaking and formation processes in the bulk triplet superfluid 
in neutron stars is investigated with taking into account of anomalous weak interactions. We 
consider the problem in the BCS approximation discarding Fermi-liquid effects. In this approach we 
derive self-consistent equations for anomalous vector and axial- vector vertices of weak interactions 
taking into account the 3 i-2 — 3 -^2 mixing. Further we simplify the problem and consider the 
pure 3 P2 pairing with rrij = 0, as is adopted in the minimal cooling paradigm. As was expected 
because of current conservation we have obtained a large suppression of the neutrino emissivity in 
the vector channel. More exactly, the neutrino emission through the vector channel vanishes in the 
nonrelativistic limit Vf = 0. The axial channel is also found to be moderately suppressed. The 
total neutrino emissivity is suppressed by a factor of 1.9 x 10 _1 relative to original estimates using 
bare weak vertices. 



1 



I. INTRODUCTION 



Thermal excitations in superfluid baryon matter of neutron stars, in the form of broken 
Cooper pairs, can recombine into the condensate by emitting neutrino pairs via neutral weak 
currents It is generally accepted that, for temperatures near the associated superfluid 
critical temperatures, emission from pair breaking and formation (PBF) processes dominates 
the neutrino emissivities in many cases. Recently [2], it has been found however that the 
existing theory of PBF processes based on the bare weak vertices violates conservation 
of vector weak current. Correct evaluations including anomalous interactions has shown 
the neutrino emission by a nonrelativistic singlet superfluid is substantially suppressed. 
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Consistent estimates of the inhibition factor can be found in Refs. 
of neutrino emissivity from the 1 Sq PBF processes was studied also in Refs. 
these are controversial (see discussion in Refs. 

Quenching of the neutrino emission found in the case of 1 5'o pairing leads to higher tem- 
peratures that can be reached in the crust of an accreting neutron star. This allows to explain 
the observed data of superbursts triggering [lO] which was in dramatic discrepancy with 
the previous theory of the crust cooling. Numerical simulations of the neutron star cooling 
in the minimal scenario 11] have shown that the suppression of the PBF processes in the 



crust of a neutron star has a significant effect at early times (t < 1000 years) and results in 
warmer crusts and increased crust relaxation times. 

We now turn to the PBF neutrino emission from the bulk superfluid neutron matter which 
is mostly caused by the triplet neutron pairin g. N eutrino energy losses due to the triplet PBF 



processes have been initially derived in Ref. 



121 ] . ignoring the anomalous weak interactions. 



From analogy with the singlet case it is clear that conservation of the vector weak current is 
violated in this approach and thus the neutrino emission in the vector channel, as obtained 



in Ref. 



is a subject of inconsistency 



131 ] . Moreover, in the triplet superfluid, the 



order parameter is sensitive also to the axial weak field. Therefore the self-consistent axial 
response of the triplet superfluid must incorporate the anomalous contributions in the same 
degree of approximation as the vector response. This effect is not investigated up to now. 

In present paper, we perform the corresponding self-consistent calculation. Formally, our 
approach is a development of Larkin-Migdal-Leggett theory [ijj], 2] to the triplet case. 
However, we discard residual particle-hole interactions because the Landau parameters are 
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unknown for a dense asymmetric baryon matter. Another reason is that the influence of the 
particle-hole interactions is not very significant in the PBF processes 4]. 

The paper is organized as follows. The next section contains some preliminary notes. We 
discuss the order parameter and the quasiparticle propagators for the triplet pair-correlated 
system with strong interactions. We also recast the standard gap equation to the form 
convenient for consideration of the processes occuring in a vicinity of the Fermi surface. In 
Sec. Ill, we formulate the set of BCS equations for calculation of the anomalous vertices and 
correlation functions of the triplet superfluid Fermi liquid at finite temperature involving a 
mixing of the 3 P2 and 3 i<2 channels lf|, ljj]. In Sec. IV, we present the general expression 
for the emissivity of the neutron PBF processes formulated in terms of the imaginary part of 
the current-current correlator. The widely used expression for the neutrino emissivity caused 
by the triplet pairing of neutrons was obtained in Ref. 12] with the aid of the Fermi golden 



rule. Therefore before proceeding to the self-consistent calculation of the neutrino energy 
losses, in Sec. V, we reproduce this formula using the calculation technique developed in our 
paper so that an apposite comparison with Ref. 12] can be made. In Sec. VI, we consider 
the anomalous vertices and the self-consistent superfluid response both in the vector and 
axial channels. Here we focus on the 3 P2 pairing with rrij =0, as is adopted in the minimal 
cooling paradigm ll|. Finally, in Sec. VII, we evaluate the self-consistent neutrino energy 
losses from the PBF processes in the triplet neutron superfluid. Section VIII contains a 
short summary of our findings and the conclusion. 

In this work we use the standard model of weak interactions, the system of units h = c = 1, 
and the Boltzmann constant ks = 1. 



II. PRELIMINARY NOTES AND NOTATION 
A. The order parameter and Green functions. 

The order parameter, D = D a/3 , arising due to triplet pairing of quasiparticles, represents 
a 2 x 2 symmetric matrix in spin space, (a,/3=t, ■!)• The spin-orbit interaction among 
quasiparticles is known to dominate in the nucleon matter of a high density. Therefore it is 
conventional to represent the triplet order parameter of the system D = J2i m ^im^a™ 
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as a superposition of standard spin-angle functions of the total angular momentum (j,rrij), 
$ S Wj) ( n ) = (2 2 a ^ Sms ) ( slm smi\jmj) Y l>mi (n) . (1) 

m s +mi=mj ^ ' 

For our calculations it will be more convenient to use vector notation which involves a 
set of mutually orthogonal complex vectors b [mj (n) defined as 

b/ mj (n) = -^Tr (g*$ jlmj ) , (2) 

where a = (<7i, <3"2, 03) are Pauli spin matrices, and g = ia^- We will use the normalization 
condition 



J 47T J 



Zmj — Su'S mjm '.. (3) 



If the most attractive channel of interactions is assumed in the states with s — l,j — 
2, 1 — j ± 1 (in the case of tensor forces) the order parameter can be written in the form 

D (n) = A ^ (d"b, mj ) g . (4) 

Irrij 

We are mostly interested in the values of quasiparticle momenta p near the Fermi surface 
p ~ pf, where the partial gap amplitudes, Ai m . (p) ~ Ai mj (pp) , are almost constants, and 
the angular dependence of the order parameter is represented by the unit vector n = p/p 
which defines the polar angles (8, tp) on the Fermi surface. 

The ground state @ occurring in neutron matter has a relatively simple structure (uni- 



tary triplet) [16j, |l7j |: 



^2A lmj b lmj (n) = A b (n) , (5) 

Irrij 

where A is a complex constant (on the Fermi surface), and b (n) is a real vector which we 
normalize by the condition 

£l» = l. (6) 
Thus the triplet order parameter can be written as 

D (n) = Ab<r£ . (7) 

We will use the adopted graphical notation for the ordinary and anomalous propagators, 
as shown in Fig. 1. 
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FIG. 1: Diagrams depicting the ordinary and anomalous propagators of a quasiparticle. 



The analytic form of the propagators can be found in the standard way 18], 19| , using 
the general form ([7]) of the gap matrix. Since the matter is assumed in thermal equilibrium 
at some temperature, we employ the Matsubara calculation technique. Then 

G {p m , p) = aG (p m , p) 5 a p , G~ (p m , p) = aG~ (p m , p) S aj3 , 

(p m , p) = aF (p m , p) hag , F {2) (p m , p) = aF (p m , p) g&h , (8) 

where a ~ 1 is the usual Green's-function renormalization constant; p m = iir (2m + 1) T with 
m = 0, ±1, ±2... is the Matsubara's fermion frequency, and the scalar Green's functions are 
of the form 

G( Pm ,p) = 2 p2 , G (p m , P ) = - , 
fm ' p rtn ' p 

F(p m ,p)= ■ (9) 
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Here 



Pi PF / x nn x 

£p= 2M^"2ilF"M^ (P ~ PF) ' (10) 

with M* = Pf/Vf being the effective mass of a quasiparticle. The quasiparticle energy is 
given by 
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E p = \lel + -TrD (n) £>t („) = ^ + A 2 6 2 , (11) 



where the (temperature-dependent) energy gap, Ab (n), is anisotropic. Here the fact is used 
that, in the absence of external fields, the gap amplitude A is real. 

Green functions of a quas iparticle (jSJ) involve the renormalization factor a ~ 1 indepen- 



dent of Co>,q, T (see e.g. [l9j). The final outcomes are independent of this factor therefore 
we will drop the renormalization factor in order to shorten the equations by assuming that 
all the necessary physical values are properly renormalized. 

The following notation will be used below. We designate as L x ,x (w, q; p) the analytical 
continuation onto the upper-half plane of complex variable u of the following Matsubara 
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sums: 



L X X' (^n,P+^;P-|) = T ^2 X (Pm + ^n,P+!) X' (p m ,P~) , (12) 

m 

where X, X' e G, F, G~, and u n = 2mTn with n = 0, ±1, ±2.... 

It is convenient to divide the integration over the momentum space into integration over 
the solid angle and over the energy according to 

d 3 p [dnl f°° , 

and operate with integrals over the quasiparticle energy: 

Zx X , (w, n, q; T) = - J de p L xx , [u,p+-,p--J . (14) 

These are functions of oo, q and the direction of the quasiparticle momentum p = pn. Here 
and below p = p F M*/ir 2 is the density of states near the Fermi surface. 

The loop integrals (|14p possess the following properties which can be verified by a straight- 
forward calculation: 

1-G-G — ZgG~ j 1-GF = —T-FG , -^G~F = -^FG~ j (15) 

1 G - F + X FG = ^Iff , (16) 

1-G-F ~ -E-FG — ~ -^^FF ■ (17) 

For arbitrary u, q, T one can obtain also 

X GG - + b 2 X FF = A+ U " ~j* Y)2 Z F F , (18) 

where v is a vector with the magnitude of the Fermi velocity V F and the direction of n, and 

A (n) ee [l G -a (n) + b 2 (n) 1 FF (n)] _ q=0 . (19) 

In the case of triplet superfluid the key role in the response theory belongs to the loop 
integrals X FF and (X GG ± b 2 X FF ) . For further usage we indicate the properties of thise 
functions in the case of u > and q — > 0. A straightforward calculation yields 

f°° ds 1 E 

and 

(X GG + b 2 l FF ) q ^ = 0, (21) 
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(1 GG - b 2 X FF ) q ^ = -2b 2 l FF (to, 0) . (22) 
The imaginary part of X FF arises from the poles of the integrand in Eq. (TSUI) at oj — ±2E: 



lml FF (u > 0, q = 0) = (cu 2 - Ab 2 A 2 ) , ^ A tanh ^- . (23) 

uy/u 2 -4PA 2 4T 

where (x) is Heaviside step function. 



B. Gap equation 

The block of the interaction diagrams irreducible in the channel of two quasiparticles, 
r a 0,7<5 5 is usually generated by the expansion over spin-angle functions (JT]). Using the vector 
notation, the most attractive channel of pairing interactions with j = 2 can be written as 



pT a ^ s (p, P') = ~ E Vlv & ( b ' m , ( n )^) (^K'm, ( n ')) > (24) 

where V/// (p,p') are the corresponding interaction amplitudes, and \l — l'\ < 2 in the case 
of tensor forces. 

In vector notation the set of equations for the triplet partial amplitudes Ai mj is of the 
form 



(P) = - E ^ / d ^ 2 ^' (P' P') A W / t" b ^ ( n >( n/ ) T E • (25) 



where 

b(n) = i^4,^(n), (26) 



A 

Inn 



as defined in Eq. ([5]). These equations can be reduced to the standard form [17| with the 
aid of the identity 

TV - 1 = ^^-tanh^M, (27) 

^?4 + £ 2 , 2£(p') 2T { ! 

and the relation 

^Tr ($ jlmj & jVmj ) = b lmj (n) ■ b?,^ (n) . (28) 

We are interested in the processes occuring in a vicinity of the Fermi surface. Therefore 
we now recast the gap equation to the more convenient form. We notice that 

' G (p m , p) G~ (p m , p) + b 2 F (p m , p) F (p m , p) , (29) 



V 2 + E 2 

fm 1 p 



i.e. Eq. ( )25|) can be written as 

V» = -^E / ^V 2 ^(p, P ')A( P ') J ^K, mj (n')b(n') 

x T E i G P') G ~ (P^» P') + PF (P"" P) F P')] • ( 3 °) 

m 

To get rid of the integration over the regions far from the Fermi surface we renormalize 
the interaction as suggested in Ref. [if]]: we define 

VP (p,p' ; T) = V W (p, j/) - V W (p,p') (GG-) n V{p (p,p'; T) , (31) 

where the loop (GG~) n is evaluated in the normal (nonsuperfluid) state. In terms of 
the gap equation becomes 

T^2[GG--(GG-) n + PFF] . (32) 



x 

J Pm,P' 

m 



and we may everywhere substitute V,j, for Vw provided that at the same time we understand 
by GG~ element the subtracted quantity GG~ — (GG~) n \(GG~) n is to be evaluated for 
oj = 0, q = in all cases]. From now we will do this and drop the superscript r on Vffi. 

Since the function GG~ + b 2 FF decreases rapidly along with a distance from the Fermi 
surface, we may replace Eq. (132]) with 

^ = -~X V ^ I ^Kv,(n)b(n)i [ dpp ^ [GG~ + b*FF] pmtp , (33) 

assuming that in the narrow vicinity of the Fermi surface the smooth functions 
Aimj (p) , Vw {p,p') , A (p') may be replaced with constants: A (p) ~ A (p F ) = A, ect.. 
The function ffT9l) is now to be understood as 



A (n) -> [l G -a - X (G - G)n + b 2 l FF ]^ =0 , (34) 
and the gap equations (|33l become: 

A Zm . = -A £ Vfo / ^b^ (n)b(n) A (n) . (35) 

The function ( J34l) can be found explicitly after performing the Matsubara's summation: 



, / n 1 f i I 1 i Ve 2 + A 2 6 2 1,6. 

A(n) = - del -. tanh tanh — . (36) 

1 ^ 4 7-oc I v^ 2 + A 2 6 2 2T £ 2T ' 
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III. EFFECTIVE VERTICES AND THE CORRELATION FUNCTIONS 



The field interaction with a superfluid should be described with the aid of four effective 
three-point vertices shown in Fig. 2. 



■ iii 

/^\ 



FIG. 2: Diagrams of the ordinary and anomalous vertices for the quasiparticle interacting with the 
external field shown by the dash line. 

There are two ordinary effective vertices corresponding to creation of a particle and a 
hole by the field that differ by direction of fermion lines. We denote these 2x2 matrices 
as f q) = r Q( g(n;o;,q) and f~ (n;cu, q) = Tp a (— n;u, q), respectively. The anomalous 
vertices correspond to creation of two particles or two holes. We denote these matrices as 
(n;u;, q) and (n;u, q), respectively. 



Given by the sum of the ladder-type diagrams 14(, the anomalous vertices are to satisfy 



the Dyson's equations symbolically depicted by the graphs in Fig. 3. 



■ i i i i 
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FIG. 3: Dyson's equations for the anomalous vertices. The shaded rectangle represents the pairing 
interaction. 

Analytically the equations reduce to the following (we omit for brevity the dependence 



of functions on u and q): 

O n ) = E(* b <->)£UE^ 



lirii 



I [ lGG -° ( &h Q f{1) - lFF { &h Q gf^ (fib) 



-X GF (<xb) (d-bj/^J f+I FG - (&K m .) ( &h ) [9T 9) 

O n ) = E($* b W n )LE v « 

-Tr 



(37) 



a/3 ■ 

dn' 



x / — it [l G - G (&b Vm .) gfM - X FF (&b Vm .) (6-b) T^g (ab) 

+X G - F (ab Vm .) gr'g (ab) - X FG (&b Vmj ) (^b) f] n , . (38) 

To obtain these equations we used the identity gg = — 1 and a cyclic permutation of the 
matrices under the trace signs. 

In general, the ordinary effective vertex is to be also found by ideal summation of the 
ladder diagrams incorporating residual particle-hole interactions. Unfortunately, the Landau 
parameters for these interactions in asymmetric nuclear matter are unknown therefore we 
simply neglect the particle-hole interactions and consider the pair correlation function in the 
BCS approximation. Thus, if the 2x2 matrix in spin space £ (n,/c) is some vertex of a free 
particle, the ordinary vertices of a quasiparticle and a hole in the BCS approximation are 
to be taken as: 

f (n,A;) = i (n,k) , f- (n,k) = £ T (-n,*) . (39) 

Discarding the particle-hole interactions, we nevertheless assume that the "bare" vertices 
are properly renormalized [ljj] in order to get rid of the integration over regions far from 
the Fermi surface, e 2 v ^> A 2 . As mentioned above, we omit the renormalization factor 
everywhere. 



G' = 




FIG. 4: Correction to the ordinary propagator of a quasiparticle in external field. 

Variation of the Green function of a quasiparticle under the action of external field U 

5G 



G' 



5U' 



(40) 
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is given by the graphs [19| shown in Fig. 4, and can be written analytically as 

G' = GG f + FF (&b)gf-g(&b) 

+ GF f {1) g (<xb) + FG (<xb) gf (2) , (41) 

where GG = G (p m + u n , p + q/2) G (p m , p - q/2), ect. 

The medium response onto external field is given by the pair correlation function which 
can be found as the analytic continuation of the following Matsubara sum 

n-K,g) = T^/g T r(fG') . (42) 

m 

IV. GENERAL APPROACH TO NEUTRINO ENERGY LOSSES 

The PBF processes are kinematically allowed thanks to the existence of a superfluid 
energy gap, which admits the quasiparticle transitions with time-like momentum transfer 
k = (w,q), as required by the final neutrino pair: k = k\ + k 2 . We consider the standard 
model of weak interactions. After integration over the phase space of escaping neutrinos 
and antineutrinos the total energy which is emitted into neutrino pairs per unit volume and 



time is given by the following formula (see details, e.g., in Ref. 20] ): 

where M v = 3 is the number of neutrino flavors; Gf is the Fermi coupling constant; and 
(x) is the Heaviside step-function. n^ ak is the retarded weak polarization tensor of the 
medium. 

In general, the weak polarization tensor of the medium is a sum of the vector- vector, axial- 
axial, and mixed terms. The mixed axial-vector polarization has to be an antisymmetric 
tensor, and its contraction in Eq. fj43j) with the symmetric tensor k^k v — ^g^ vanishes. 
Thus only the pure- vector and pure-axial polarizations should be taken into account. We 
then obtain Im n^ ak ~ Cy Im Ely' + C\ Im n^, where Cy and C*a are vector and axial- vector 
weak coupling constants of a neutron, respectively. 

V. PRESENT STATUS OF THE PROBLEM 

The widely used expression for the neutrino emissivity caused by the triplet pairing of 



neutrons was obtained in Ref. 



12J with the aid of the Fermi "golden" rule. Therefore before 
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proceeding to the self-consistent calculation of the neutrino energy losses, it is instructive 
to reproduce this formula using the calculation technique developed in our paper. We will 



prove the result of Ref. [12| can be obtained from our equations ( 1431) and (1421) if to remove 
the field interactions through anomalous vertices [second line in Eq. (1411) ]. We will label the 
corresponding results with tilde. 

The authors of Ref. [12j state that the weak current of nonrelativistic neutrons is caused 
mostly by the temporal component of the vector current, Jo — and by the space 

components of the axial-vector current, j = \E' + o"j\I / . Consequently to reproduce their 
result we need to evaluate the temporal component of the polarization tensor in the vector 
channel and the spatial part of the axial polarization. Omitting the anomalous contributions 
for the temporal component of the vector polarization we have to substitute for 

T = T~->i, f (1 ' 2) -» , (44) 

where 1 is a unit 2x2 matrix in spin space. Eq. ( 14~2]) is valid for each of the tensor 
components. Inserting the temporal component of the vector vertex, from Eqs. f T4T]) . ( |4~2"]) 
we then obtain after a little algebra: 

n™(ui, q )=4 P J^ 1 -(l GG -b*l FF ) . (45) 

In obtaining this expression we used Eqs. (fl~2l) . ( TT4l) and the identity (<rb) (<rb) = b 2 . 

Only small transferred momenta, q < u> ~ T, contribute into the neutrino energy losses. 
Since the transferred momentum comes in the polarization function in a combination qV F <C 
u), A (Fermi velocity V F is small in a nonrelativistic system), to the lowest accuracy, we may 
evaluate the polarization tensor in the limit q = 0. (In the same approximation the above 
authors evaluate the matrix elements of a quasiparticle transition.) Then using Eqs. ( 122|) 
and ( )23l) we find 

ImfT°>>0,g = 0) = -47rp f ^p A 2 -^fc££L tanh — . (46) 

Polarization tensor in the axial channel can be evaluated in the same way. In this case, 
omitting the anomalous contributions we have to take 

f (n,Jfe) &i , f~ (n,Jfe) -> of , -> . (47) 

Then we find after some algebraic manipulations: 

flf (u, q) = Ap J ^ f~ {1 GG - b 2 X FF ) Sa + Iff (48) 
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In obtaining this we used the identities ga- T g = <x, and & (^b) = 2b — (<xb) er. 
With the aid of Eqs. ([22]) and ([23]) we find: 

Im («, > 0, t = 0) = -*r„ J - [S, - -tf J fA' ^_^J tanh - (49) 

Inserting the imaginary part of the polarization tensor into Eq. f[43]) we calculate the con- 
traction of ImII^ ak with the symmetric tensor kpk, v — k 2 g^ u to obtain 

= -4vrp / — 6 2 A 2 ) _ ! tanh — 

J ^ 2w V / w 2 -46 2 A 2 4T 

x [C 2 (g 2 + gi)+C 2 (2(a; 2 -g 2 )-gi)] , (50) 
where gy and q± are defined as 

tf-pW.al = «•-,}■ (51) 

After a little algebra we obtain the neutrino emissivity in the form: 



= 9M^ PfM * [ —A 2 r du; , - 1 (C 2 + 2C 2 ) , (52) 

120«* PF J 4vr n 7 2An (i + exp# ) 2 v / u; 2_4 A 2 ^ 

where A n = A b (n). 

With the aid of the change u = 2T^Jx 2 + A 2 /T 2 one can recast this expression to the 



form obtained in Ref. 



12|: 



_4G|A/; 2 , 7 /-rfnA 2 r z 4 

where z = ^ x 2 + A 2 /T 2 . 

Apparently the contribution of the vector channel in this expression is a subject of incon- 
sistency, since conservation of the vector current in weak interactions requires cully (u, q) = 
qjly (uj, q) , and thus one should expect Ily (to > 0, q = 0) = for the correct result instead 
of Eq. (1461) . This however was not proved explicitly for the case of triplet pairing. We now 
focus on this calculation. 
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VI. ANOMALOUS CONTRIBUTIONS 



A. Vector channel 

The self-consistent longitudinal polarization function ITy {oj > 0, q) incorporates the 
anomalous contributions. At finite transferred space momentum the problem of determin- 
ing the vertex corrections is much complicated. Typically massless Goldstone modes that 
arise due to symmetry breaking play a crucial role in conservation of the vector current. In 
the anisotropic 3 P2 phase rotational symmetry is broken and three Goldstone modes arise 



(termed angulons in Ref. 21|). However, since we are interested in the specific case of q = 
the temporal component of the anomalous vertex T M (// = 0, 1, 2, 3) can be retrieved from 
the Ward identity which requires j^|: 

wT (1,2) (n;w, q) - qT (1,2) (n;w, q) = ±2D (n) . (54) 



From this identity we immediately find 



2A- 

T w (n;u;,q=0) = — bag, (55) 



and 



T (2) (n;u;,q = 0) = - — gba . (56) 

0J 



In the BCS approximation, the ordinary scalar vertices are to be taken, as given by Eq. 
3J). Inserting the above vertices into Eqs. (jUJ, ( l4"2j) we obtain after a little algebra: 



n° (c,q = 0)=4p|^Q(X GG -6 2 X FF ) +^ 2 X TO ) • (57) 
Using Eqs. dHJ), (JTZD yielding 



we finally find 

n°° ( W> q = 0) = 2C 2 v p J ^ {lac + b 2 l FF ) q=0 . (59) 

Comparing this with Eq. ( 12 ip we obtain ITy (ou, q = 0) = 0, as is required by the current 
conservation condition. We found that the neutrino emissivity through the vector channel 
vanishes in the limit q = 0. This proves explicitly that the neutrino emissivity via the vector 
channel, as obtained in Eq. (1531) . is a subject of inconsistency. 
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B. Axial channel 



We now focus on the axial channel of the weak polarization. The order parameter in the 
triplet superfluid varies under the action of axial-vector external field. Therefore the self- 
consistent axial polarization tensor also must incorporate anomalous contributions. Then 
from Eqs. (14T|) . ( 142|) we obtain after simple algebraic manipulations: 



Iff M = 4p 



(in 

Air 



\ {Igg - b 2 l FF ) 5 l3 + VT FF h -h 
03 l FF ^Tr (oilfg (<rb) - &i (<rb) gff 



2A 



(60) 



As in above, we focus on the case q = and omit for brevity the dependence on n and u. 
The anomalous axial- vector vertices Tr' 2 ' (j = 1, 2, 3) are to be found from Eqs. ( 1371) . ( |38l) . 
where the ordinary vertices are given by Eq. ( |47lh 

Up to this point we have not discussed the n dependence of bi mj (n). This makes Eq. 
f )59|) valid in the case of tensor forces resulting in the 3 P 2 — 3 F 2 mixing, because the general 
form of Eqs. ( )37|) . ( 138|) for the anomalous vertices takes into account not only spin-orbit 
interactions but the tensor interactions in the channel of two quasiparticles. Now we simplify 
the problem according to approximation adopted in simulations of neutron star cooling 11 1 
and consider the case of paring in the 3 P 2 channel, when 1 = 1, and Vw = SwV, and the 
vectors b mj (n) are given by 

bo = ^ (-ni, -n 2 ,2n 3 ) , 

bi = -b*_i = - (n 3 , m 3 , n x + in 2 ) , 

b 2 = b* 2 = y|("i +m 2 ,mi -n 2 ,0) , (61) 

where ni = sin^cos<y9, n 2 = sin 9 sin ip, n 3 = cos 6. From now on we will drop the subscript 
I = 1 by assuming b m . = b 1>TOj , A m = A 1>m ., etc. 

We will focus on the p-wave condensation into the state 3 P 2 with rrij = which is 
conventionally considered as the preferable one in the bulk matter of neutron stars. In this 
case, Eq. (jSJ) implies 

b (n) = b (n) , A = A (62) 
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and the gap equation (I35I) reads 



1 = -Vj ^P(n)A(n). 



(63) 



From Eqs. (1371) and fl38|) we obtain the vertex equations of the following form (i = 1, 2, 3): 



da! 



I GG -Tr g *b*)T 



-I FF Tr[[&b* )(*b) gf® (6-b 



=.(1) 



-—X FF 2i[h\, I) 



(64) 



^ (2) (n) =V^§&b^(n) I ^ [z G - G Tr ((<xb mi ) g% 
- l FF Tr ((<rb m ,) (6-b) ff^ (o-b)) 



-—X FF 2i (b mj xb). 



(65) 



In obtaining the last line in these equations we used & (6-b) = 2b — (6-b) & along with 
Tr ((crb^) (6-b) 6") = 2* (b^xb), and Eqs. (USD, (USD- 

Inspection of the equations reveals that the anomalous axial- vector vertices can be found 
in the following form 

fW(n lU ) = 5]BWMK.)9, (66) 



T< 2 > (n, W ) = $>$M$(*b^) 



(67) 



These general expressions can be simplified due to the fact that the function X FF (n; a;) given 
by Eq. (|20|) is axial - symmetric, and the last (free) term, in Eqs. AM]) and (|65|) . can be 
averaged over the azimuth angle to give 



and 



dip 
2^ 

dtp 



b* x B) = / g (b; x B) 



\/6r, 



^(b%xb)=0 

2tt V - 2 ' 



(bj x b) = 



(b*i x b) = -e 



(68) 



(69) 



2tt v 1 ' 4 ' 7 2tt v _1 y 4 

where e = (1, — i, 0) is a constant complex vector in spin space. The following relations can 

be also verified by a straightforward calculation: 



b* b m /. = <L,. m /.b* b 



2tt m * m > 



mjm'j'-'mj'-'rnj i 



(70) 
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dip 
2^ 



bb 



bb m ' 



3 j 



( bb m 3 ) ( bb ^) 



Relations and allow to conclude that Bg 1,2 -* = B^ 2 '* = 0, and 

T« (n)= (b« (<xbi)+B«(<xb-i))s, 

tW(n) = y(BP(*bI)+Bg(W_ 1 )) . 

Inserting these expressions into Eqs. ( jMJ) and ( |65|) . taking the traces and using the 
onality relations (j3J) along with relations (JTOjl . ( ITT)) , and 

b 2 = b 2 , b^b! = b^b-i , 
(bbl) (bbi) = (bblj (bb_ X ) , 

we obtain the equations: 



B 



(i) 
±i 



X gg -B2 (bxK 



-Z^Bg ((b*b x ) b 2 - 2 (bjb) (bbx)) - ^F F e^P 



and 



B 



(2) 
±1 



4:71 



-X FF Bg ((b x bl)6 2 - 2 (b^) (bbl)) + ^ F e*^6 2 



Solution to Eqs. (171)1 . ( 175)) can be found in the form 

b« = -bl 1 ; , bL 2 ; = -b«, 

where 

B x = e / (w) , B_! = e* / (w) 
and the function / (co) satisfies the equation 
dn 



f = -v 



4tt 



(X GG - + b?X FF ) (bJbO / - 2X FF (bjb) (bbi) / - ^FF^-b 2 



Using Eq. (II 8p we can rewrite this as 



f = -v 



dn 

47T 



.4 



^Iff) (Kb,) f - 21 FF {bib) (bbx) / - ^ Fi ^6 2 
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At this point it is convenient to recast the left side of this equation according to Eq. (|63|) : 



f=-Vf 



In this way we obtain the equation 
(in 



/ 



47T 

3 u 
22A 



An 



{bib, - b 2 ) A + 2 ( (bjbx) - (bjb) (bbx) ) 1 FF 



-V (n A n 



(80) 



4A 



cfn T 



47T 



6 2 X, 



(81) 



Since the function Xpp (n; u) is axial - symmetric and 

& 2 = i(l+3n*) , btb, = | (1 + nl) 



(bb*) (bbx) 



-n. 



(1 " nl) 



(82) 
(83) 



Eq. flHTl) can be integrated over the azimuth angle, yielding the following solution 



/ (w,<7 = 0) 



^ vi ^ r *4 f 1 + 3n ') z - r » • (84 » 



where 



X [u, q = 0) = / drz 3 
./o 

3 / u 1 

+ 



i (l-3ng) A (n 3 ,T) 



4 V 2A 2 



(l + Tig) - fig (l - Tig) J 1 F F (™3, w, T) 



(85) 



and the functions A (rig, T) and Xp^ (rig, w, T) are given by Eqs. (1361) and (!20|) . 

Explicit evaluation of Eq. fl84|) for arbitrary values of a; and T appears to require nu- 
merical computation. However, we can get a clear idea of the behavior of this function 
using the angle-averaged energy gap in the quasiparticle energy, <( A 2 6 2 ) = A 2 . (Replacing 
angle-dependent quantities in the gap equation with their angular average has been found 
to be a good approximation In this approximation the functions Xff{^, T) and A (T), 

in Eqs. (1841) and fl85l) . can be moved beyond the integrals. Using also the fact that 



A [ dn 3 (l - 3n 2 3 ) = 
Jo 



we find 



/ 



Au 
2w 2 - A 2 /5 



(86) 



(87) 
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Thus, in approximation of the average gap, the function / (co) is real-valued and is indepen- 
dent of the temperature. 

Poles of the vertex function correspond to collective eigen-modes of the system. Therefore, 
the pole at u 2 = A 2 /5 signals the existence of collective spin oscillations. The decay of the 
collective oscillations into neutrino pairs gives the additive contribution into neutrino energy 
losses. However, examination of the collective modes deserves a separate study, which is 
beyond the scope of this paper. Here we concentrate on the PBF processes discussed in the 
introduction. 

In this case tu > 2Ab (6) > \^2A and, to obtain a simple analytic approximation, we 
omit a small term A 2 /5 in the denominator of Eq. (j87p . thus obtaining the axial- vector 
anomalous vertices in the following simple form: 



(88) 



f ( 2 ) (n) = yjl^g (e (crbO + e* (*b_0) , (89) 

Having obtained this simple result we can evaluate the axial polarization function. In- 
serting dSSJ) and dSHD into Eq. (ED]) gives 



nX(«) = 4/»/£ 



\ (x GG - l 2 x FF ) 5 l3 + l 2 x FF h -h 



+ - Si-sdjs) -b X FF 

' q=0 

The first line in Eq. ( 190]) can be evaluated with the aid of Eq. fl22|) . We find: 



(90) 



= - A pjt ~ W ~ I {5ij " W ) PlFF (w ' q = 0) ■ (91) 

Using Eq. ( 123]) we obtain the imaginary part of axial polarization: 

Im Tlf (u > 0, q = 0) 

dn / hbj 3 /r r r ,\ 72 ,6 (uj 2 — 46 2 A 2 ) w 
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VII. SELF-CONSISTENT NEUTRINO ENERGY LOSSES 



As we have obtained Imlly (cu > 0, q = 0) = 0, using Eqs. (123|) and (f9~Tj) we find 



T2x2 Q(u 2 -4b 2 A 2 ) , w 
x fe 2 A 2 — v , = 



tanh ^- . (93) 



Contraction of this tensor with {k^k v — k 2 g^ u ) gives: 

- M < 2 - * - a / ^ 2 ^fcS tanh £ • (94) 

where 

ql = l{^)\ql = q 2 -q 2 . (95) 

The rest of the calculation is already performed in Sec. VII. The neutrino energy losses 
can be written immediately after inspection of Eqs. ( 150]) and ( 194)) . From this comparison it 
is clear that in order to obtain the correct neutrino energy losses, it is necessary to replace 
the factor (C^ + 2C\) with (1/2) C\ in Eq. (EH}. In this way we obtain 



2 



G 2 F C\K PF M^ [ ^ r V * , 2 , (96) 
157T y 47r y (l+expz) 



where A 2 = A 2 5 2 (n) = |A 2 (1 + 3 cos 2 9), and z = a/x 2 + A 2 /T 2 . Comparison of this 
expression with Eq. (I55|) shows that the neutrino energy losses caused by the 3 P2 pairing in 
neutron matter are suppressed by the factor 

1 r 2 

~ 0.19 (97) 



2 (C 2 + 2(7 



with respect to that predicted in Ref. 12] 



For a practical usage we reduce Eq. (f96|) to the traditional form 

where M and M* are the effective and bare nucleon masses, respectively; c is speed of light, 
and 

f dn A 2 f 00 , 2 4 

F t = / — — £ / rfx 5 . 99 

J 4vrT2 y (1 + expz) 2 1 J 
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Notice the gap amplitude A (T) denned above is a/2 times larger than the gap amplitude 
Aq (T) used in Ref. 12| , where the same anisotropic gap A n is written in the form A 2 = 



Aq (1 + 3 cos 2 9). However, the function F t , defined in Eq. (I99j) . is independent of the 
particular choice of the gap amplitude, therefore the analytic fit (B) suggested in Eq. (34) 



of Ref. 12), is valid and can be used in practical computations. 



VIII. SUMMARY AND CONCLUSION 

In this paper we have performed a self-consistent calculation of the neutrino energy losses 
due to the pair breaking and formation processes in the triplet-correlated neutron matter 
which is generally expected to exist in the neutron star interior. Since the existing theory of 
anomalous weak interactions in the fermion superfluid is well developed only for the case of 
1 S'o pairing we have generalized the corresponding equations for the triplet pairing including 
the case when the attractive tensor coupling is operative. 

Exact solution of the vertex equations is much complicated because of anisotropy of the 
triplet order parameter. Fortunately only small values of the transferred space momenta are 
significant for the considered processes in the nonrelativistic approximation. Therefore the 
weak vertices as well as the polarization functions can be evaluated in the limit q = 0. 

Before proceeding to the self-consistent calculation we reproduced the neutrino energy 



losses as obtained in Ref. [12], using the calculation technique developed in our paper. We 
have shown that the result of Ref. 12[ can be obtained in the BCS approximation from our 
equations (l4~3l) and ( H2l) if to remove the field interactions through anomalous vertices. 

The exact solution we found for the vector part of the weak polarization, 
riy (u > 0, q = 0) = 0, is consistent with the current conservation condition. This general 
result, which is obtained including the tensor couplings and the Fermi-liquid interactions, 



means that the neutrino emissivity in the vector channel, as obtained in Ref. 12j, is a 
subject of inconsistency. 

The self-consistent consideration of the axial weak polarization is more complicated. 
In this case, inclusion of the tensor forces and the Fermi-liquid effects requires numerical 
computations even in the limit of q = 0. Therefore to obtain a simple analytic result we have 
considered the 3 P2 pairing in the state with rrij = which is conventionally considered as the 
preferable one in the minimal cooling scenario of neutron stars. We have also neglected the 



21 



residual particle-hole interactions since the Landau parameters are unknown for the neutron 
matter at high density. 

Finally we used the self-consistent polarization functions for evaluation of the neu- 
trino energy losses due to PBF processes in the 3 P 2 neutron superfluid with rrij = 0. 
The obtained self-consistent neutrino emissivity, is given by Eq. (l9T)j) . This expression 
needs to be compared to the emissivity (153|) originally derived in Ref. 12], ignoring the 
anomalous weak interactions. One can see the neutrino emissivity is strongly suppressed 
due to the collective effects we have considered in this paper. The suppression factor is 
(1/2) Cl/(C$ + 2Cl)~ 0.19. 

Since the neutron 3 P 2 pairing occurs in the core, which contains more than 90% of 
the neutron star volume, the found quenching of the neutrino energy losses from the PBF 
processes can affect the minimal cooling paradigm. 
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